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abstract 

We investigate real-time tachyon dynamics of unstable D-brane carrying fundamental string 
charge. We construct the boundary state relevant for rolling of modulated tachyon with gauge 
fields excited on the world-volume, and study spatial distribution of the fundamental string 
charge and current as the D-brane decays. We find that, in contrast to homogeneous tachyon 
rolling, spatial modulation of the tachyon field triggers density wave of strings when electric 
field is turned on, and of string anti-string pairs when magnetic field is turned on. We show that 
the energy density and the fundamental string charge density are locked together, and evolve 
into a localized delta-function array (instead of evolving into a string fluid) until a critical time 
set by rolling tachyon's initial condition. When the gauge fields approach the critical limit, the 
fundamental strings produced become BPS-like. We also study the dynamics via effective field 
theory, and find agreement. 
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1 Introduction 



Recently, there has been considerable progress in understanding decay of unstable D-brane in 
string theory The decay proceeds, as analyzed first by Sen [1], via rolling of the tachyon 
living on the D-brane world-volume from the top of the tachyon potential to the bottom (closed 
string vacuum). It was found that, in the weak coupling limit g st = 0, as the tachyon rolls down 
the potential hill, the unstable D-brane is converted to a pressureless gas, referred as "tachyon 
matter". Details of the decay process is of considerable physical interest, as it would shed light 
on semi-classical and non-perturbative aspects of string dynamics. 

An interesting physics question is, as the D-brane decays, how various constituents bound 
inside an unstable D-brane are liberated. To answer the question, we will consider a bound- 
state of unstable D-brane with fundamental string, and study dynamical evolution of energy and 
charge distributions. In the previous works [10, 15], this problem was studied for homogeneous 
rolling of the tachyon, and it was found that the constituents form a uniformly distributed 
fundamental string fluid on the hyper-surface of the D-brane world- volume, behaving BPS-like 
in the limit we take the world-volume gauge fields to a critical value. In this work, we shall 
examine whether and, if it does, how spatial modulation of the rolling tachyon field would affect 
fate of constituent's distribution in the ambient space-time. In approaching the problem, we 
shall continue utilizing the approach set forth in the previous work [15]. 

We begin with salient features concerning homogeneous rolling of the open string tachyon. 
The real-time dynamics is described [1], in the boundary state approach, by turning on an 
appropriate boundary interaction to the c = 1 conformal field theory of the X°-field: 



AS = X J dt j> da 8(t) cosh X°(t, a) 



Boundary state corresponding to the above interaction was constructed in [1], and takes the 
form 



|D25) T(x o) = \B) x o \N) xi ® (ghost), 



where 



\B) x o = /(£>>+ ^V-i^-ilO) 



+ 
+ 



hix^a^-cf, + V-i) 2 + ((a°_ 1 ) 2 a°_ 2 + a^ct,) 2 )] |0> 



\N) xi = exp -£-«-n7SLn |0> 



71=1 



n 



(* = !,•••, 25) 



;i. 2 ) 



t Among various works contributed to the progress, a list of those more relevant to the present work in- 
cludes [1]-[18]. 



I ghost) = exp \ - (b-nC-n + 6- n C- n )J (Co + Co)ciCi|0) . 

The coefficient functions f(x°), g(x°), hi(x°), h 2 (x°), h,3(x°), ■ ■ ■ describe time evolution of the 
rolling tachyon in the closed string channel, and are computable explicitly, as was done in 
[1, 10, 11]. Expanding the boundary-state Eq.(l.l) in oscillator levels, one can obtain explicit 
form of the linear coupling of the tachyon-rolling D25-brane to each closed string mode. It was 
found that, at late time, the rolling tachyon evolves into a pressureless "tachyon matter". It was 
further found that the coupling to higher-level closed string states is hierarchically pronounced 
[11], a feature escalated even further in a situation the world- volume electric and magnetic 
fields are turned on [10, 15]. It also entailed an issue on the stability of the "tachyon matter" 
once a non-zero g st is turned on. 

The rolling of spatially modulated tachyon field can also be studied using the boundary 
state approach [7]. In [7], Sen considered the inhomogeneous tachyon corresponding to the 
following boundary interaction: 

da cosh —= cos — =. 

It was found that half of the total energy start to accumulate at localized points and form an 
array of delta-function-like singularity at finite time. Furthermore, after the system hits the 
singularity, the total energy suddenly reduces to zero and the energy density vanishes eventually 
everywhere. 

In this paper, we turn on gauge fields along with the inhomogeneous rolling tachyon, and 
study time evolution of the decaying D-brane. Surprisingly, we find that, in the presence 
of spatial modulation of the tachyon field, magnetic field (as well as electric field) induces 
fundamental string constituents on the D-brane world- volume. We also find that time evolution 
develops a singularity not only for the energy density but also for the fundamental string charge 
density. Consequently, our study indicates the following picture of the D-brane decay. Initially, 
driven by the tachyon field modulation, the fundamental string energy and charge densities 
evolve into free streaming of thin flux-tubes, whose thickness approaches zero up to a critical 
time set by the initial condition of the tachyon and gauge fields. The critical time is the moment 
a singularity develops to the boundary state — in particular, energy and charge densities are 
discontinuous across the critical time, leading allegedly to energy and charge non-conservation. 
Physically, as the energy and the string charge ought to be conserved, a possibility is that 
macroscopic fundamental strings of zero thickness are built up until the critical time and then 
are decoupled from the system right after the critical time (recall that we are working with g st 
set to zero.). 

We study couplings to the massless as well as some of the lower-level massive closed string 
modes, which can be extracted from the boundary state. Through these couplings, the bound- 



ary state acts as a space-time varying source for the closed string field in the linearized equation 
of motion 



where Qb is the BRST operator and | <& ) is the closed string field. We find that the coupling 
to the closed strings generally blows up sharply as the system evolves across the critical time, 
suggesting that the source | B) of the closed string field equation in Eq.(1.3) may be approxi- 
mated as a sort of impulse. Then, once g st is turned on, the equations of motion for both open 
and closed strings will be altered by nonlinear interactions, and our analysis may cease to hold 
close to the critical time. The situation is analogous to the formation of a heavy star in general 
relativity: once the energy density becomes huge compared to the Planck scale, one cannot 
ignore the nonlinear terms in the Einstein field equation. 

Several previous works [20, 21, 22] examined fate of the fundamental string constituents 
using the effective field theory approximation to an unstable D-brane, and claimed that a 
diffusive string fluid of arbitrary density profile is formed at the tachyon potential minimum, at 
least at the classical level. Our boundary state (as well as effective field theory) analysis shows 
the contrary that thin string-like flux tubes can be formed via real-time tachyon rolling *. 

This paper is organized as follows. In section 2, we begin with recapitulating Sen's construc- 
tion of the boundary state representing rolling of modulated tachyon in bosonic string theory, 
and examine couplings to massless and massive closed string modes. We present detailed anal- 
ysis of the time evolution of these couplings, and indicate that the couplings to both massless 
and massive closed string modes exhibit a singularity at the critical time. In section 3, adopting 
the recipe prescribed in [15], we construct the boundary state describing rolling of modulated 
tachyon with constant electric and magnetic fields. We extract time evolution of energy and 
charge density associated with fundamental string and analyze the behavior. In section 4, we 
study the tachyon rolling dynamics again, but now in the effective field theory approach. We 
find similar but somewhat differing result compared to those obtained from the boundary state 
approach. We contrast our results to earlier works [20, 22], and assert that stability criterion 
used in these works needs to be reconsidered in our case. Section 5 is devoted to discussion 
including some speculation concerning the physics around the singularity. 

2 Rolling of Modulated Tachyon 
2.1 Boundary State Construction 

Following [7] , consider the boundary-state description of the rolling of inhomogeneous tachyon 
on a D25-brane in the bosonic string theory. In the boundary state approach, it amounts to 
tThis might be related to the result of [13]. 
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turning on the following boundary action: 

/X° X 1 
da cosh —= cos —= (2.4) 
V 2 v 2 

to the c = 2 conformal field theory (CFT) associated with X^X 1 of the bosonic D25-brane 
[23, 24]. The corresponding boundary state describes a D25-brane whose world- volume tachyon 
field is modulated sinusoidally in space x 1 , and rolling in time x°. Schematically, it takes the 
form: 

|D25) T(x o iX i ) = |S) x o )X i(g)^ 2 |7V) Xi (g,|ghost). (2.5) 

The boundary action Eq.(2.4) turns out to define an exactly solvable perturbation, as the c = 2 
CFT boundary state \B) x o t x 1 is obtainable in a factorized form: 

\B) x0 , x i = \B)+®\B). , (2.6) 

where \B)± are the c = 1 boundary state Eq.(1.2) with the replacement of X° — > (X°±iX 1 )/\/2 
and A — > A/2 [7]. Explicitly, we find that 

\B) ± = fH^^'M 

+ ^g ± (x°,x 1 )(a°_ 1 ±ia 1 _ 1 )(cP_ 1 iialJlO) 

+ i/if ir 1 )^ 2 ± ^i 2 )(«° 2 ± iai 2 )|0) 

+ S 1 )^ x ± iaii) 2 ^! ± ^!) 2 |0) 

8 

+ ^(s ,* 1 ) [(a^±iaij 2 (z£ 2 ±iai^ |0) 
+ ••• 

where the coefficient functions are defined by 

^(x^x 1 ) = M +e + ^^sin(A7r/2)J + M + e VT~ s in(A7r/2) J -1, (2.7) 
g ± (x°,x 1 ) = 2cos 2 C\n/2)-f ± (x°,x l ), (2.8) 
h±(x°, x 1 ) = 2 cos 2 (Att/2) - 2 sin(A7r/2) cos 2 (Att/2) cosh ( ^-^- ^J - /±(a; , x 1 ), 

^(x ,* 1 ) = 4sin(A7r/2)cos 2 (A7r/2)cosh( j + ^{x , x 1 ), 

hfix^x 1 ) = -2sin(A7r/2)cos 2 (A7r/2)sinh (^-^^J . (2.9) 

Note that the functions labeled with + superscript are complex conjugate of those labeled 
with — superscript. Note also that, except hf(x°,x 1 ), all coefficient functions are related to 
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Expanding the matter part of the boundary state Eq.(2.5) in powers of the oscillators and 
coupling functions associated with them, we obtain 



|D25) T = F(s°,a;)|0) 



^Iabcdix^x 1 ) a^a^aiiolilO) + % - J ab c(x°, x l ) (a a _ 1 a b _ 1 a c _ 2 + a^a^at^ |0) 

(2.10) 



+ 4 



where the coupling functions are given by 



F=\\f + \\ 2 , 

Goo = -Gn = Re(/V), Goi = +Gi = M/V), 

^ = -||/ + ||% (for i,j = 2,3,-..,25), (2.11) 

for the tachyon and massless level modes, and 

#oo = -H u = Re(/+/0, #io = +#oi = Im(/ + /ir), 

Hij = -\\f + \\% (for i,j = 2,3,---,25), 

/oooo = /mi = +Re(/ + /i2) + ||s + || 2 , 

/oon = /noo = -Re(/ + /i2 ) + ||# + || 2 , 

/oioi = +/1010 — +/0110 — +/1001 — —R e (f + h 2 ), 

/oiOO — ~ /oill = +/0001 — — /llOl = +I m (/ + ^2 )) 

/ijfej = ||/ + || 2 (^fe^ + SuSjk) (for i, j, /c, Z = 2, 3, • - • , 25), 

/001 = -/in = Re(/ + /i 3 ), Jooo = -/no = -Im(/ + /i3 ), 

/010 = +/100 = Mf + K ), Jon = +J101 = +Im(/ + /i3 ), (2.12) 

etc. for the first massive level modes. Various relations among the coupling functions are 
attributable to the specific factorized form of the boundary state, as given in Eq.(2.6). 

2.2 Time Evolution of Coefficient Functions 

As alluded above, the D25-brane boundary state | D25 ) T acts as a source to the closed string 
field equation Eq.(1.3). The source being space-time varying, we thus need to examine its time 
evolution. In Eq.(2.10), the source is decomposed into closed string mass- levels, and we shall 
be interested in the evolution of the coupling functions Eqs.(2.11, 2.12). 

As f ± (x° ) x 1 ) are the building blocks for most of these coefficients, we will analyze their 
behavior first. From Eq.(2.7), we see that they are singular at x°±ix 1 = =F\/2 log (— sin(A7r/2)) . 



To be specific, consider the case < sin(A7r/2) < 1 in the following. Flipping the sign of A, 
interpreted as flipping the sign of the initial value of the tachyon field, amounts to relocating 
the x 1 -locus of the singularities by \/2n. When < sin(A7r/2) < 1, the singularities are located 
at x 1 = y/2(2n + l)n (n G Z) at the critical time x° = ±x®, defined by 



x° c = v^log ( 1 ) . 

H|sin(A7r/2)|y 



(2.13) 



Let us examine closely the behavior of f ± (x°,x 1 ) around the critical time. Useful formulae 
for our analysis are: 

lim ( l -— + -=) = l±2V2nY 5 (x - V2(2n + l)n) , (2.14) 



and 



/■2\/27r 

Jo 



dx 1 



1 _ for |a| > 1 

1 + ae ±ixyV2 ~ \ for a\ < 1 



(2.15) 



The behavior of the real part of f ± (x°,x 1 ) close to the critical time x° = x° c is obtainable 
by making use of Eq.(2.14): 

= ±2y/2n ]T 5 (x l - y/2(2n + l)n) 
nez 

+ [(l + sin 2 (A7r/2)e i:cl/ ^) _1 + (l + sin 2 (A7r/2) e - i:cl/ ^) _1 - l| . (2.16) 

Note that the sign of the 5-function part (the first line) flips as x° passes through the critical 
time x° c , whereas the regular part (the second line) remains the same. Note also that, for a 
fixed x° ~ x° c , the real part Eq.(2.16) is an even function of x 1 across each singularity located 
at x° = \[2{2n + The imaginary part of / ± (a; ,x 1 ) is also singular at the critical time x° c : 



Jm (/ + (xV)-/~(sV 

■ ( x 1 \ 
—i tan — = 

\2V2J 

(l + sin 2 (A7r/2)e" 1 /^) _1 - (l + sin^A^e^/Vy 1 



(2.17) 



In this case, there is no discontinuity in x° for both the singular part (the first line) and the 
regular part (the second line). For a fixed x° ~ x° c) the imaginary part Eq.(2.17) is an odd 
function of x l across each singularity at x° = \[2{2n + l)ir. 
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0.25 



Figure 1: Behavior of the energy density (be Re/ + (x°, x 1 )) (vertical axis in arbitrary unit) for 
X = 0.3 in x = [0, 2] and x 1 = [0, in the left, and its section just before the critical time 



x c ~ 1.116 ■ ■ ■ in the right. 



It is also of interest to examine the spatial averages of the functions f ± (x°, x 1 ). We obtain 
them from Eq.(2.15) as 



<rx* u ) := 



1 



dx 1 f^x ,! 1 ) 



for \x I > x c 

1 for < x c r . 



(2.18) 



2\/27r Jo 

Generically, functions associated with massive modes, such as hf,hf, depend on f ± (see 
Eq.(2.9)), and hence develop the same singularity as / at the critical time x° c . There are, 
however, some exceptions and some other functions, such as h 3 , do not depend on and 
evolve in time in a regular manner. 



2.3 Coupling to Closed String Modes 

Having obtained the coupling functions Eqs.(2.11, 2.12), let us now examine time evolution of 
each of them. 

• graviton coupling 

The energy-momentum tensor of the decaying D25-brane can be read from the boundary state 
Eq.(2.10) as [1] 

^(x ,^ 1 ) = ^ (Gabix ^ 1 ) - riabFix ^ 1 ] 
Utilizing the expressions Eq.(2.11), we obtain [7] 



T o = -T n = T 25 cos 2 (A7r/2)Re/ + (a; ,x 1 ) 
Toi = +T 10 = T 25 cos 2 (Avr/2)Im/ + (x°,x 1 ) 



-T 2 ,\\f+{x\x l )\\ 2 
7 



(for i = 2, 3, •■•,25). 



(2.19) 



Figure 2: Behavior of the momentum density (be Im/ + (x°, x l ) ) (vertical axis in arbitrary unit) 
for A = 0.3 in x° = [0, 2] and x 1 = [0, 2\^2n] in the left, and its section at the critical time 
x° c ~ 1.116 ■ ■ ■ in the right. 

One readily sees that the energy-momentum tensor behaves as a sort of impulse source to the 
closed string field equation. From Eq.(2.18), the average energy is 

1 f 2 ^ , n m , n is f for |ac°| > x° c 

o • 

c 



2^71" ./o 1 T 25 cos 2 (A7r/2) for < x 



It indicates that the energy suddenly disappears after the critical time! This comes about as 
follows. As can be seen from Eq.(2.16) and Eq.(2.18), half of the total energy is squeezed 
at x 1 = y / 2(2n + 1)tx (n E Z) and form the 5-function singularities as the critical time x° c is 
approached, while the other half is smoothly spread out along the x 1 -direction (See Fig.l). After 
the critical time x° c , the sign of the 5-function flips (see Eq.(2.16)) and cancels out contribution 
of the other half that are smoothly spread out. The momentum density is zero on the average, 
but develops a singular profile which may be approximated by ~ tan(a; 1 /2v^2) (See Fig.2). 

To gain a better picture concerning the origin of the singularity, let us examine the energy- 
momentum conservation carefully. Introduce for convenience complex coordinates z = (x° + 
ix 1 ), z = (x° — ix 1 ) so that do = (d z + b\), d\ = i(d z — c\). The coefficient functions / are 
analytic functions of z,~z, respectively. From Eq.(2.19), we then find that 

8 a T aQ = -r 25 cos 2 CXn/2)[d z f-(z) + d I f + (z) 

d a T al = zr 25C os 2 (Xir/2)[d z f-(z)-c\f + (z)\. (2.20) 

As alluded in the previous section, f + {z), f~{z) have simple poles at z = — \pl log(— sin(A7r/2)) 
and z = +v^2 log(— sin(A7r/2)), respectively. Thus, the first relation in Eq.(2.20) is proportional 
to Snez 3( x ° ~ x < ^.)5(x 1 — v / 2(2n + l)vr), invalidating the energy-momentum conservation. On 
the other hand, the second relation in Eq.(2.20) vanishes, yielding the correct force law. 
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Figure 3: Behavior of closed string tachyon and dilaton coupling density (be \\f + (x°,x 1 )\\) 
(vertical axis in arbitrary unit) for A = 0.3 in x° = [0, 2] and x l = [0, 2\/2%]. 

Physical interpretation of such singular behavior is not clear to us, so we shall primarily 
concentrate our consideration on the early evolution before hitting the critical time, i.e. \x°\ < 
x° c . As the total energy ought to be conserved, a plausible possibility suggested in [7] is that 
the missing energy has escaped through the singularity to form an array of codimension-one 
D-branes. We leave viability of this scenario aside for future study. (See section 5 for further 
discussion.) 

• tachyon coupling 

The tachyon coupling density is given by 

Ptach y on(x°,X 1 )=r 25 ||/ + (x ,X 1 )|| 2 , 

and its behavior is plotted in Fig. 3. The coupling again develops a singularity at the critical 
time x° = x° c . From Eqs.(2.11,2.12), we also find that couplings to all higher closed string 
modes with polarization along i — 2, 3, - • • , 25 directions behave the same way as the tachyon 
coupling. 

• dilaton coupling 

As the part of the boundary state that couples to the sigma model dilaton is 

(c + c ) (c_iCi + C_ X Ci) |0) gh , 
the dilaton coupling of the decaying D25-brane is given by 

Pdi i,Ux°,x 1 )=T 25 \\f + (x ,x 1 )\\ 2 . 

Its behavior is the same as that of the tachyon coupling and the pressure T# (i — 2, 3, • • • , 25). 

• massive string mode coupling 

Coupling to massive closed string modes, as compared to tachyon and massless string modes, 

9 



Figure 4: Behavior of —R,e(f + h 1 ) in the left and lm(f + h 1 ) in the right (vertical axis in 
arbitrary unit) with A = 0.3 in x = [0,4] andx 1 = [0,2y/2n]. 

entails certain new features. As is seen from Eq.(2.12), these couplings are distinguished from 
those for tachyon and massless modes in that, in addition to f ± (x°, x 1 ), new coefficient functions 
hf (x°, x 1 ), hf (x°, x 1 ), hf (x°, x 1 ), ■ •• are involved. We thus expect that time evolution of the 
massive mode couplings would display qualitatively significant departure from those of the 
tachyon and the massless state couplings. Generally, they do not vanish at late time, though 
the energy density vanishes soon after the system hits the singularity at the critical time. As an 
example, let us consider H a t,(x°, x 1 ) couplings. They are governed by the real and the imaginary 
parts of the coefficient functions f ± and hf: 

H 00 + tH 01 = -H n + tH 10 = f + K . 

We plot the real and the imaginary parts of the function 

f+h^ = 2cos 2 (Att/2) 

in Fig. 4. The plots show that a singularity develops precisely at the critical time x° = 
where the couplings to tachyon and massless modes diverged. This is again attributable to the 
fact that all these couplings depend on the function f ± . On the other hand, for the asymptotic 
behavior, the H a b coupling differs from the couplings to tachyon and massless modes. At 
x° ^> x®, we find that 

f + h^ ~ -2cos 2 (A7r/2)e- v/ 2 ia;1 , 

exhibiting persistent modulation along the x 1 -direction. Note that, as compared to the early 
epoch x° the modulation wave- number has increased twice asymptotically. 

The coupling functions at higher massive levels would be governed by two qualitatively 
different contributions. These coupling functions are extractable from expanding the product 
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2 sin(A7r/2) cosh 



[x 



IX 



V2 



f + -\\r 



of two c = 1 boundary states Eq.(2.6) in powers of closed string oscillator levels. For the first 
massive level modes, we have seen in Eq.(2.12) that the couplings are quadratic products among 
the coefficient functions (f ± , g ± , hf, hf, hf) up to level-2. Likewise, for a generic massive level 
mode, the coupling would be quadratic products among the coefficient functions, now involving 
higher-level ones. Expanding the products, we see that the coupling function comprises of two 
kind of contributions: the one depending on (e.g. most of Eq.(2.12) in the case of the first 
massive level modes) and the one independent of f :± (e.g. products among the first parts of 
g ± ,hf,hf, hf in the case of the first massive level modes). As the first type of contributions 
depends on f ± , there will always be a singularity at x° = x° c . The second type of contributions 
is non-singular, generally increasing as i° — > oo, and hence is qualitatively analogous to the 
behavior of the couplings for homogeneous tachyon rolling. 

A remark is in order. So far, we have focused primarily on constructing the boundary state 
Eq.(2.10) itself. It might well be that some components of the boundary state Eq.(2.10) are 
BRST exact, viz. couple only to off-shell or unphysical closed string states. In that case, 
recalling that the boundary state Eq.(2.10) acts as a source to the closed string field equation 
Eq.(1.3), the BRST exact components would not be contributing to on-shell closed string 
processes, such as closed string radiation out of decaying D-brane. Whether there indeed 
exists BRST exact components and, if so, which modes (and their time-dependent coefficient 
functions) belong to them are pertinent questions to be examined in case one needs to study 
on-shell processes. In this work, we will leave the issue aside for future study. 

3 Turing on Electric and Magnetic Fields 

With the physics motivation addressed in section 1, we extend the analysis of section 2 and 
construct a boundary state for an unstable D-brane with electric and magnetic fields turned 
on. 

3.1 Rolling of Modulated Tachyon under Electric Flux 

We shall now consider turning on a constant electric field on the D25-brane world-volume, 
and study rolling of the spatially modulated tachyon field. As is well-known, the electric field 
induces fundamental string constituents on the D25-brane world-volume. If the tachyon field 
were spatially homogeneous, the induced string charge density (which is given by the electric 
displacement field) would be homogeneous over the D25-brane world-volume as well. On the 
other hand, if the tachyon field were spatially modulated, the induced string charge density 
would be modulated accordingly and are driven to evolve as the tachyon rolls down to the 
potential minimum. In such a circumstance, an interesting and potentially important question 
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is whether the displacement field flux can possibly get squeezed into (an array of) thin flux tubes. 
If so, adding spatial modulation on a decaying D-brane would be a new way of manufacturing 
a macroscopically large, fundamental closed string! 

With such motivation, consider turning on the world- volume electric field along 2-direction 
e = Fq2. With tachyon modulation along 1-direction, this amounts to a fluid of fundamen- 
tal strings, whose macroscopic configuration is stretched along 2-direction and charge density 
varies across 1-direction. As prescribed in the work [15], the corresponding boundary state 
is obtainable by a sequence of T-dual map, boost by e, and inverse T-dual map, all along 
the 2-direction. Denote the new coordinates as y a, s and the new string oscillators as /3°'s. 
Starting from Eq.(2.10), the new boundary state would again be expressible as: 



\B) 



F e (y°y)\0) 



+ G^(yV)/£i/£i|0> + ^(vV)/£2/£ 2 |0> 



+ 
+ 



^(vV)/£i/£i/£i/£i|o> 
\K* c {y\ y 1 ) {P-1P-1P-2 + P-2P-1P-1) |o> + 



where 



oitivW) = r'^-'G^jrV.y 1 ) 

H e ab (y°y) = 7- 1 (*A- 1 ^A) ak ( 7 -y,y 1 ) 
ZtaivW) = t-'Ca-^'A-^A) A) ( 7 - 1 y°,y 1 ) 



etc., and 



A 



7 7 e 
1 

7e 7 



V 



7 



J 



Explicitly, the new coefficient functions are given as 



Ggo(l/V) 
G?i(i/V) 



r'FirV.y 1 ) 

7- 1 ll/ + (7-V,y 1 )H 2 , 

7[Goo(7'V,y 1 )-e 2 G 22 (7-V,y 1 ); 
-r l \\f + {T l y\y l )\\ 2 + 2 7 cos 2 (A7r/2)Re/ + ( 7 -V,y 1 ), 

7- 1 Cii( 7 -y,y 1 ) 
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= +7' 1 ||/ + (7-V,y 1 )l| 2 -2 7 - 1 cos 2 (A7r/2)Re/ + (7-V,l/ 1 ), 

GUy\y l ) = 7[G 22 (7"V,l/ 1 )-e 2 Goo(7-V,l/ 1 ); 

= -7- 1 ||/ + (7-V,y 1 )l| 2 -27e 2 cos 2 (A 7 r/2)Re/ + (7-V,y 1 ) ) 

= -7- 1 ||/ + (7-V,y 1 )H 2 (for < = 3,...,25), 

G e 02 (y ,y 1 ) = -G e 20 (y°,y 1 ) = 7 e [G 00 (i~ V, y 1 ) ~ G 22 ( 7 ~ V, y l ) 

= 2 7 ecos 2 (A7r/2)Re/+(7-y,?/ 1 ), 

C^(yV) = +G?o(l/V) = Ci^W) = 2 cos 2 (A7r/2) Im/ + (7 _1 y°, y 1 ), 

GiafoV) = -G e 2l {y\y l ) = eG^vW) = 2ecos 2 CXn/2)lmf + ( 1 - 1 y°,y 1 ) 

for the tachyon and the massless level modes, and 



7 



H, 



oo 



e 2 H 22 



"7 



H e 0Q (y°y) 

^TiG/V) = r l H u = -T l Ke{f + h- 1 ){r 1 y\y 1 ) 
HUy°y) = 7- 1 ^« = -7- 1 ll/ + (7- 1 y°,y 1 )ll 2 

H^(y ,y 1 ) = +H! (y°,y 1 ) = H 01 = Imtf+Wfr-V, V 1 ), 

H e 12 (y°,y 1 ) = -H e 2 i(y°,y 1 ) = eH 01 = eim(/ + / i r)(7"V,y 1 ), 

Re(/ + /0 + ll/ 



He(/ + /ir) + e 2 ||/ + || 2 (7-Y.y 1 ), 



e 2 Re(/ + / i r) + ||/ + || 2 l (7- 1 Z/°,y 1 ), 



(for z = 3, •• -,25), 



H e 02 (y ,y 1 ) = -H e 20 (y°,y 1 ) = ei [H 00 - H 22 ] 
for the coupling to the first massive level modes. 



e7 



(7~W) 



3.2 Coupling to Closed String Modes 

Coupling to the closed string modes are readily obtained from the boundary state, as we did 
in section 2. 

• tachyon and dilaton couplings 
For tachyon and dilaton fields, the couplings are given by 

pU yo «(y°y) = p^U*/ ^ 1 ) = t 257 - 1 ||/ + (7-V,i/ 1 )II 2 • 



We see that, as the electric field is turned on along x 2 -direction, the rolling is time-dilated and 
the spatial modulation along ^-direction remains unaffected. In the critical limit, e — * 1, these 
couplings vanish because the overall Born-Infeld factor 7" 1 dilutes the coupling density. 

• graviton coupling 
The energy-momentum tensor is given in terms of the coupling functions by 

T^ivW) = Tj f (Gl ab) (y o y)- Vab F%y°y)) , 
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so we obtain 



+ T 257 cos 2 (A7r/2)Re/ + (7-yV), 
+T 25 cos 2 (A7r/2) Im/ + ( 7 - V, y 1 ), 
-T^- 1 cos 2 (Att/2) Re/ + ( 7 - 1 t/°, y 1 ), 

-r 25 e 2 7 cos 2 (A7r/2)Re/ + (7-V,y 1 )-r 257 - 1 ||/ + ( 7 -V,y 1 )|| 2 , 
-Tvn^Wf+irVrfW (for < = 3,...,25). (3.21) 

Let us compare Eq.(3.21) with the energy-momentum tensor Eq.(2.19) at zero electric flux. 
Apart from the expected Born-Infeld time dilation and Lorentz contraction, we see that the 
pressure along the electric field direction receives a new contribution, proportional to the energy 
density. This part is attributed to the fundamental string constituents. 

• Kalb-Ramond coupling 
At nonzero electric flux, coupling to the Kalb-Ramond field 

<&(M ■= f^O/V), 

is newly induced, whose non-zero components are 

QpfllG/V) = T 257 ecos 2 (A7r/2)Re/ + ( 7 - 1 2/ °, 2/ 1 ), 

Q[i2](yV) = T 25 ecos 2 (A7r/2)Im/ + ( 7 - 1 2/ ,y 1 ). (3.22) 

As anticipated, the induced Kalb-Ramond coupling is proportional to the nonzero electric field, 
thus making up the decaying D25-brane to carry a fundamental string charge. Note that the 
charge density of the fundamental string stretched along 2-direction is given by Q[o2], and it is 
proportional to the energy density in Eq.(3.21) as 

Q m (y ,y 1 ) = eT* (y°,y 1 ) . (3.23) 

Likewise, the current density of the fundamental string Q[ 12 ] is proportional to the momentum 
density in Eq.(3.21) as 

Q [ i2 ] (y°,y 1 ) = eT* 1 (y°,y 1 ) . (3.24) 

Therefore, just as the behavior of energy density explained in section 2.3, half of the fundamental 
string charge density would squeeze to form (an array of) localized 5-function profile as the time 
y° approaches the critical time y° c . The off-diagonal (Ol)-component of the energy-momentum 
tensor and (12)-component of the Kalb-Ramond coupling imply that the fundamental string 
density is actually flowing along the 1-direction! This is to be contrasted with the situation 
of homogeneous tachyon rolling, where turning on electric field e induces a fundamental string 



T 2 e 2 (2/V) = 
^(2/ V) = 



14 



fluid on the decaying D-brane world- volume, but the fluid is at rest. As shown in [15, 16], 
turning on magnetic field b perpendicular to the electric field lets the fluid to flow rigidly. We 
now find that modulated tachyon field can also trigger the fundamental string gas to flow. 
Moreover, the flow velocity varies spatially, albeit it is determined by the modulation of the 
rolling tachyon itself. 
• BPS Limit 

Now that the electric field e is restricted to |e| < 1, Eqs.(3.23, 3.24) imply the following BPS-like 
inequalities 

Tq > \ Qfo2]\ and \Tq ± \ > |Q[ 12 ]| , 

which are saturated precisely at the limit |e| — > 1. 

One important effect of turning on the electric field is the time dilation by the 7 _1 factor: 
x° — > 7 _1 y°. Therefore, the critical time is now dilated to 

y° c =iV2hg( , 1 ) = 7 xg. 

V| sm(A7r/2)|/ 

If we take the extremal limit 

e — > 1 and A — > 1 

while holding the average energy (E) = ^T 25 cos 2 (A7r/2) fixed, we find that the critical time is 
situated at a finite value: 

E 



y° c 



V2T 25 ' 

indicating that the relaxation time scale grows with the total energy one starts with. 
In the extremal limit, the real part of f ± becomes 

Ref ± h -i o iv f 2V2nJ2nez5(y 1 -V2(2n+l)n) for \y°\ < y° c , 
I for \y°\ >y° c . 

Thus, prior to reaching the critical time y®, non- vanishing components of the energy- momentum 
tensor and the Kalb-Ramond tensor current density become 

Too = ~ T 22 = I Q02I = 2V2nE ]T 5 (y 1 - V2(2n + 1)tt) . (3.25) 

neZ 

This is the same relation as that saturated by an array of BPS fundamental strings stretched 
along the |/ 2 -direction, except that we are now considering homogeneous distribution of the 
string fluid along the transverse y 3 , ■ ■ ■ , y 25 -directions. It thus indicates that the fundamental 
string constituents are confined to co-dimension one hypersurfaces located at y 1 = y/2(2n+l)n. 
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Note, however, that this extremal limit is not a smooth limit for the full boundary state. 
As pointed out in [10, 15] for the homogeneous case, the coefficients in the boundary state for 
higher-level massive modes with more than three oscillators for 0- or 2-directions are weighted 
with extra 7 factors, and hence will diverge in the extremal limit. In light of the remark at 
the end of section 2, a possible way out is that (part of) the non-smooth modes belong to the 
BRST exact class. 

3.3 The Effect of Magnetic Field 

Physics-wise, it is also of interest to turn on magnetic field on the decaying D-brane world- 
volume. The world-volume gauge field is combined with pull-back of the Kalb-Ramond field 
X*B a b into a gauge-invariant combination (F a j ) +X*B a } ) ). It thus implies that nonzero magnetic 
field induces current density of the fundamental string on the decaying D-brane world-volume. 
Here, we examine how the fundamental string constituents are affected by the magnetic field. 
We also consider the case with pure magnetic field. In this case, the total fundamental string 
charge is zero, but, nevertheless we find that polarized fundamental string charge density is 
induced, and the system evolves to a configuration with an array of fundamental string - anti- 
fundamental string pairs (i.e. an array of fundamental strings of alternating orientations). We 
will see that, again, spatial modulation of the rolling tachyon field plays a prominent role. 
Following the prescription given in the previous work [15], we find that the boundary state 

is 

\B) T , e+b = F e+6 (yV)|o> 

+ <^(yV)/£i£iio> + 

+ ••• , 

where the coefficient functions are obtained as 

Git'iy ^ 1 ) = r'r 1 ('(fiA)- 1 GAfi) afc ( 7 "V,ry -^A 

Kt\y\y l ) = 7- 1 r 1 C(^A)- 1 ^Afi) a6 (7'V,rV-te7i/ ), 

C c b d (y°y) = 7- 1 r 1 C( fiA )'H'^ A ) _1/Afi ), c Afi ) arf ^~V,rV-fe7i/ ), 
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etc., in terms of transformation matrices 

\ 



A = 



f 7 7e' 

10 

7e' 7 

V 



and 



n = 



/ 1 o o 

7 —76 
+76 7 



and 



7 



1 



,/2 



with e' 



and 



7 



We then find that the nonzero components of the energy-momentum tensor are 



rpe+b 
J 00 

rpe+b 

1 m 



rpe+b 
1 11 



+T 25 77(1 + b 2 ) cos 2 (Att/2) Re/ + , 
+T 25 cos 2 (Att/2) flm/ + + 77e6Re/" 

-T 25 cos 2 (A7r/2) 



77(l-e 2 )Re/ + -- 



2e6 



Im/~ 



T 2 e + b = +T 25 cos 2 (A7r/2) 



+ 6 2 ) 

77(& 2 -e 2 )Re/ + + IT ^ 2T Imf 
2 (for i = 3, • - • , 25), 



while the nonzero components of the Kalb-Ramond coupling tensor are 

)e+b 

'[02] 



Qf n +f = T 25 cos 2 (A7r/2) [77c Re/ + - 6 Imf 



O e+6 

^[12] 



T 25 cos 2 (Att/2) 



e[ TT ^j Im/ + + 77&Re/- 



where Z^ 1 := / ± (7 _1 y°, ^~ l y l +ebqy°). One can check that the conservation equations d a T^ b = 
and d a Q e g fr b = are satisfied. Couplings to other closed string modes are also obtainable 
analogously from the prescription given in [15]. Before proceeding further, we will contrast the 
above results against those obtained in the previous sections. 

First, the spatial coordinate x 1 is now replaced by (7" 1 ?/ 1 + eb^yy ) in the functions f ± , so 
the whole system would be moving in the ^/-direction with a constant velocity 



V := 



eb 



< 1 . 



1 + b 2 

Indeed, prior to the critical time the average momentum 

1 /•2v / 2tt „ 

(P) := —=- / dy^iyW) = 77e6T 25 cos 2 (A7r/2) 
2V27T J 



(3.26) 
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and the average energy 

(E) := / dy 1 T£ b (y°, y 1 ) = 77(1 + b 2 ) T 25 cos 2 (A7r/2) 

2v2vr Jo 

are nonzero, and the velocity of the energy flow is computed as (P) / (E) = eb/(l + b 2 ). This 
is consistent with Eq.(3.26). 

Second, we note that the (01)-, (11)-, (22)-components of the energy-momentum tensor are 
now composed of both real and imaginary parts of f ± . Recall that real and imaginary parts 
of / ± are even and odd functions, respectively, of the comoving spatial coordinate (y 1 + Vy°), 
and develop a singularity at the critical time This gives rise to several effects. For (01)- 
component of the energy-momentum tensor, Poynting field momentum of the electromagnetic 
fields (proportional to eb in magnitude) adds up into an inhomogeneous distribution along in- 
direction, and accumulates to an array of delta functions as the critical time y® is approached. 
For the (11)- and (22)-components of the energy-momentum tensor, the electromagnetic fields 
induce additional contribution proportional to lmf + whose magnitude is determined by the 
velocity V. 

Third, [02]- and [12]-components of the Kalb-Ramond coupling tensor indicate that the 
fundamental string charge and current densities are induced not only by the electric field e 
but also separately by the magnetic field b. We are thus led to a surprising observation that, 
even without the electric field, a decaying D25-brane can carry fundamental string constituents 
induced by the modulated tachyon and the magnetic flux! In this real and imaginary 

parts of f^ 1 are even and odd functions of (y 1 + Vy°), the net charge would be zero while the 
net current is nonzero — the system is a sort of plasma of the fundamental strings (consisting 
of strings of alternating orientations). 

Of particular interest is the extremal limit. Take 

e -> Vl + b 2 and A — > 1, 

while keeping the average energy (E) held finite. Then, the terms which survive in this limit 
are 

T 2 (2y/2nE6(yW)) = — T 2 T 2 e 2 +6 (y°, y 1 ) 
T 2 v(2V2nE5(y°,y 1 )) = V T^ b (y° , y 1 ) 
TV 2 (2V2nES(y°y)) = V'T^y^y 1 ) 
-2V2itE5(y\y l ) 
2V27iET5(y°,y 1 ) 

2V27rEVT5(y°,y 1 ) = VQ?+\ (y°, y 1 ) . 



n +b (y°,y l ) = 

GKf(vV) = 
Q&KvW) = 
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Here, V = b/Vl + b 2 is the velocity Eq.(3.26) in the extremal limit, T = 1/y/l - V 2 = V / TTF 
is the corresponding gamma factor, E = (E) /T 2 is the average energy in the rest frame, and 

6(y°, y 1 ) := £ * + Vy°) - V2(2n + 1)tt) 

neZ 

is an array of the delta functions. As such, each component of the energy-momentum and 
the fundamental string current tensors have a constant magnitude. Evidently, if we observe 
the system in the rest frame, the results are in agreement with the BPS equality, Eq.(3.25). 
Therefore, the system is viewed most transparently as an array of BPS fundamental strings 
stretched along the y 2 -direction, but now moving in ^-direction uniformly with a constant 
velocity V. 

It would also be of interest to consider the case with large magnetic field to learn the 
effect of pure magnetic field. If we set e = and take b — > oo, the leading terms for the 
energy-momentum tensor and the Kalb-Ramond coupling tensor are 

75,(1/°, i/ 1 ) * +T 25 cos 2 (A7r/2)6Re/ + ( 1 / , 2 / 1 ) 

7*(yV) ~ +T 25 cos 2 (A7r/2)Im/ + ( 2/ °,y 1 ) 

7*(yV) * Oib' 1 ) 

7*(yV) - +T 25 cos 2 (\*/2)bRef + (y ,y 1 )-T 25 b\\f+\\(y°,y 1 ) 

T+ivW) * -T 25 &||/ + ||(</V) 
Qfo^y ^ 1 ) ~ -T 25 cos 2 (A7r/2)6Im/ + (y°,y 1 ) 

Q\ l2] (y Q ,y l ) ~ +T 25 cos 2 (A7r/2)Re/ + (y°, 2/ 1 ), (3.27) 

with / ± (|/ ,i/ 1 ) = f ± (y°,by 1 ). Now, the period in the ^-direction is shortened to 2y/2n/b. 
The energy and momentum densities are proportional to Re/ + and Im/ + , respectively, which 
is qualitatively similar to the previous results Eq.(2.19) or Eq.(3.21). On the other hand, 
the fundamental string charge and current densities are now proportional to Im/ + and Re/ + , 
respectively, which is opposite to what we have found for the case with pure electric field 
Eq.(3.22). The behavior of the fundamental string charge and current densities can be read 
from Fig. 2 and Fig.l, respectively. Since Im/ + = at y° = 0, the charge density is initially 
zero everywhere, but then starts to get polarized as the tachyon rolls down the potential hill. 

Intuitively, such a behavior can be understood as follows. As one can readily see from the 
gauge invariance X*B — > X*B + dA 1; F — > F — dA 1; the magnetic field on a D-brane world- 
volume induces fundamental string current density. However, when the tachyon field is at the 
bottom of its potential and the D-brane disappears, there is nowhere for the string current to 
flow. Therefore, the region of the D-brane world-volume where the tachyon is near the bottom 
of the potential acts just like an insulator for an electric current. In our situation, as the 
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modulated tachyon evolves, the regions around y 1 = 2^/2nir/b (n G Z) behave as insulators 
while those near y 1 = \/2(2n + l)n/b (n £ Z) behave as conductors for the string current. 
Thus, our system is analogous to a series of capacitors injected with a non-zero electric current, 
in which electric charge starts to accumulate at the edges of the conductors. The polarization 
of the fundamental string charge density is induced in this way. When the system reaches the 
critical time, the current density develops a delta-function singularity, leading to a configuration 
with a polarized fundamental string - anti-fundamental string pair near the critical time. 

4 Effective Field Theory Approach 

It is also of interest to compare the results obtained from the boundary-state approach with 
those from other approaches. In this section, we analyze the D-brane decay in the effective 
field theory approach, extending previous results [15] to the case with inhomogeneous tachyon 
rolling. Consider the world-volume Lagrangian of an unstable Dp-brane: 



£dbi = ~T P V(T)J- det( V + F) F(z), (4.28) 



where V(T) is a tachyon potential and T{z) is a function of 

/ -i \ (oft) 

d a Td b T, 



normalized such that .F(O) = 1. The energy- momentum tensor T ab and the Kalb-Ramond 
current tensor Q ab are obtainable by replacing 77^ — > i] a b + h a b, F a b — * F a b + b a b in the Lagrangian 
Eq.(4.28) and expanding in powers of h a b,b a b. 

SCdbi = |v(T)fd^+fj V ab (h {ab) + b [ab] ) + --- 

= ±(T^h (ab) + Q [ab] b [ab] )+---, (4.29) 

where 

v,t = T(z) {^fT + 2F{Z) {^hT dMT fcbf ■ (430) 

Here, we consider rolling tachyon field with spatial modulation in x 1 -direction, 

d T = f, d{T = T', d 2 T = ■ ■ ■ = d p T = 0, (4.31) 

and, for simplicity, turn off all other components of the gauge fields but F 02 = e. Then, 
Eq.(4.29) and Eq.(4.31) imply that nonzero components of the energy-momentum tensor T ab 
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are 



T 00 = +T pl V(T)(f(z)+2 1 2 f 2 r(z)), 

T m = +T pl V(T)2fT'F'(z), 

T n = -T pl - 1 V(T)(F(z)-2T' 2 F'(z)), 

T22 = -T p e 2 1 V{T)(F{z)+2 1 2 t 2 F\z))-T pl - l V{T)F{z), 

T u = -T p7 -V(T)^), (for i = 3,4,...,p), (4.32) 

and nonzero components of the Kalb-Ramond current tensor Q[ a b] are 

Q02 = eT pl V(T)(F(z) + 2 1 2 t 2 F\z)), 

Q12 = eT pl V(T)2TT'F(z), (4.33) 

with z = (-7 2 T 2 + T' 2 ) and 7 = 

In the effective field theory approach to unstable D-brane based on the Lagrangian Eq.(4.28), 
second and higher derivatives of the tachyon and gauge potential are neglected. As such, in 
case the tachyon evolves sharply and rapidly, we would not expect a good agreement with 
the corresponding results of the previous section, where the boundary state approach surely 
indicated rapid and sharp evolution of the tachyon field. Despite such shortcomings anticipated, 
we see in the above result a behavior observed in the boundary-state approach in Eq.(3.21) and 
Eq.(3.22). For example, relations 

Qo2 = eT 00 , Qi 2 = eT 01 , T 22 = -e 2 T 00 + T H . (4.34) 

featured by Eqs.(3.21, 3.22) also follow from Eqs.(4.32, 4.33). In particular, we again obtain 
from the first equation in Eq.(4.34) that the energy density is proportional to the fundamental 
string charge density. 

One might think from the above observation Eq.(4.34) that, in so far as the electric field e 
is kept fixed, any configuration is classically stable since energetics does not tell which one is 
favored among those with a fixed total fundamental string charge. However, this is not correct. 
As we have seen explicitly in the previous section, the system does evolve to thin electric 
flux tubes without changing total energy and fundamental string charge during the conversion 
process. The reason is as follows. In the conventional stability analysis, one compares the energy 
of static configurations with a fixed total charge. If one finds a minimum energy configuration 
among them, this configuration is considered stable. This is because the kinetic energy is 
always positive-definite and the minimum energy configuration cannot evolve dynamically to 
other ones without increasing the total energy. Note that, in this argument, one has implicitly 
assumed that variation of the fields does not affect the total charge. In our situation, however, 
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it is not enough to consider the energy of static configurations with a fixed fundamental string 
charge, as the kinetic energy § is also seen to contribute to the fundamental string charge. 
Even if one finds a configuration of a minimum energy among the static ones with a fixed 
charge, it is always possible for the configuration to evolve to some other configuration with 
"less" fundamental string charge, (here we mean "less" provided we set T = by hand), as 
the rest of the fundamental string charge can be carried off by the kinetic energy. This is to 
be contrasted against the reasoning of [20, 22], where the conventional stability analysis was 
taken for granted. 

If we allow the electric field e to vary, the minimum energy configuration with a fixed total 
fundamental string charge is realized in the extremal limit |e| — > 1. Note that, in this limit, 
one also needs to set V(T) — > from the outset to keep the total energy finite ^. This is a 
special limit in which the time evolution is frozen because of the (by now well-understood) 
time-dilation effect of the electric field, and can be easily seen from the fact that this system 
is T-dual to a codimension-one unstable D-brane moving with the speed of light [15]. In this 
limit, it is claimed that, since the system does not evolve at all, one can engineer a string fluid 
of an arbitrary energy density distribution [20] . On the other hand, if one starts from a generic 
point of the potential hill and consider the dynamical decay of the unstable D-brane, there is no 
reason that the system would evolve to a string fluid of arbitrary energy distribution advocated 
in [20] . Indeed, we shall now show that the system evolves to the contrary. 

Suppose we start with an inhomogeneous tachyon configuration with V(T) ^ and |e| < 1, 
as considered above. The pressure along the a; 1 -direction, pi = Tn, is, in general, non-zero and 
gets modulated along the x 1 -direction. Then, the momentum density along the a^-direction, 
T 01 , cannot remain constant, so the energy (and hence fundamental string charge density) 
distribution would change every moment during the evolution. In order to see whether the 
flux will be squeezed and form an arbitrarily thin string-like flux tube or diffuse to a uniform 
distribution, it is imperative to follow the dynamics of the system carefully. 

As an example, let us examine the background-independent string field theory action for 
non-BPS D-branes in superstring theory. We have [25] 

V(T) = e-^ 2 and T{z) = Z -^^. (4.35) 

If we take a linear-gradient profile of the tachyon field T = ux 1 and T = as an initial tachyon 
configuration, the energy density p and the pressure p along the x 1 -direction are 

p = T 00 = +T pl V{ux 1 )F{u 2 ), 
pi = T 11 = -T pl - 1 V{ux 1 )V{u% 



§more generally, contribution of time-dependence of the tachyon to the total energy 

^Here, we have assumed that the kinetic energy is positive-definite, i.e. T{z) + 2j 2 T 2 f(z) > .F(O) = 1. 
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where V(z) := T{z) — 2zT\z). In this case, we see that the ^-coordinate dependence is 
carried solely by the tachyon potential, V{ux x ). One can show that the function T[u 2 ) is 
always positive for any value of u 2 > 0, so the energy density p is manifestly positive-definite. 
The energy density p is accumulated around the origin, x 1 = 0, with a typical width of order 
0(u~ l ). It was shown in [25] that the total energy decreases monotonically with u. The energy 
density becomes delta-function distribution pi^x 1 ) oc S^x 1 ) in the u — > oo limit, representing a 
BPS D(jo— l)-brane with the electric field. The function V(u 2 ) is also a positive function with 
lim^oo V(u 2 ) = 0. Therefore, in so far as u 2 is finite, the pressure p\ is negative-definite with 
the minimum located at the maximum of the energy density. From the continuity equation 
<9 T i + dipi = 0, we see that the momentum density T 01 is turned on, and the energy starts to 
accumulate to the center. This clearly suggest that the the system will evolve to a thin string 
configuration. In order to learn detailed time evolution of the system, one will need to solve 
the equations of motion and follow the evolution explicitly. 

In the above analysis, we considered the evolution under the initial condition of a kink-like 
tachyon profile. It carries a D(p- l)-brane charge, and hence the decay product would be 
a bound state of the D(p — l)-brane and fundamental string constituents. In order to obtain 
purely fundamental string constituents in the final state, we will need to consider a configuration 
with T" 7^ 0. In this case, it is difficult to analyze the system based on the effective action 
Eq.(4.28), as the second and higher derivatives of the tachyon field are assumed small and 
thus truncated. Qualitatively, however, until the tachyon modulation has not grown up large 
enough, the energy density and the pressure of the system are well approximated by Eq.(4.32). 
Recalling that the qualitative behavior is governed primarily by the tachyon potential V, we 
expect that the pressure p reaches its minima at locations where the energy density reaches 
its maxima. As such, again, we expect that the inhomogeneity tend to grow as above. (See 
[8, 9, 17, 18] for related arguments.) 

5 Discussion 

In this paper, to gain deeper understanding concerning dynamical fate of the fundamental 
string constituents, we have studied rolling of the inhomogeneous tachyon in the presence of 
constant gauge fields on the world-volume of the decaying D-brane. Quite interestingly, we 
have found that the fundamental string charge density is induced not only by the electric field 
but also by the magnetic field. Although the gauge fields turned on are uniform, distribution 
of the fundamental string constituents is spatially inhomogeneous because of modulation of 
the rolling tachyon field. We have shown that the tachyon rolling defines an exactly marginal 
boundary perturbation to the boundary state and have solved the system exactly. It represents 
an exact classical solution of string theory, at least prior to reaching the critical time, where 
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the system develops a singularity. 

The most pressing point for further study is a physical interpretation of the critical time 
and singularity therein. Notice that our consideration corresponds to the weak coupling per- 
turbation theory at the lowest order, and we ignored radiation of massless and massive closed 
string modes, triggered by the rolling tachyon, as well as back-reaction to the decaying D-brane 
and background geometry. Indeed, we have seen that the total energy is conserved and stored 
in the open string sector until the critical time. The singularity seems to suggest that the weak 
coupling limit is not a smooth limit to take. Once we turn on g st , no matter how tiny it is, the 
coupling to closed string modes blow up at the critical time. So, the physics with g st = and 
that with tiny g st could be completely different at the critical time 

In the more realistic set-up with small but non-zero g st , inferring from Eq.(1.3), it is plausible 
to suppose that the energy will be transferred to a collection of closed strings near the critical 
time. After some energy is converted to closed string modes, it would appear as if the energy is 
lost for the observer looking only at the open string modes. It is not clear whether the energy 
is radiated away or stored at the singularity to form a macroscopic fundamental string or a 
D-brane. Once we turn on non-zero g st , since the coupling to closed string becomes impulsive 
and huge near the singularity, there ought to be a large back reaction. In particular, it is not 
hard to imagine that the space-time around the energy core will be severely curved to form 
an event horizon. If it is the case, the energy cannot escape through it, at least classically. 
This could be a possible scenario for the formation of a macroscopic fundamental string or a 
D-brane. If the string coupling g st is nonzero (no matter how small), nonlinear coupling of the 
tachyon to other open string and closed string states becomes important near the singularity. 
It would be certainly of interest to take these string states into consideration, and obtain a 
more complete picture of the rolling tachyon dynamics. 

In [19], it was claimed that the non-perturbative confinement mechanism is responsible for 
the formation of string-like object in unstable D-brane system. It was further argued in [20, 21] 
that classical dynamics has no mechanism to make the electric flux tube thin. In [22], classical 
stability analysis for the unstable D-brane effective action was made, and it was claimed that 
the flux tube is unstable against uniform spread-out. However, our result suggest that the 
flux can squeeze dynamically to a thin flux tube even at the classical level, in contrast to the 
conclusion drawn in [20, 21, 22]. Note that [20, 21] considered the situation that the tachyon 
potential is set to zero from the outset (instead of scaling limit), so it is delicate to make a 
direct comparison of their results with ours in which we considered the tachyon rolling down 
from a generic point of the potential hill. 

Furthermore, in the works [20, 22], (in)stability of the flux tube was analyzed by looking 

II Sec [11] for related issues for the late time behavior of the homogeneous rolling tachyon. 
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for a static configuration minimizing the total energy among those with a fixed fundamental 
string charge. With the identification of the correct electric displacement field, one readily sees 
that all configurations with the same fixed charge have the same total energy so it appears 
superficially that any of these configurations is classically stable. However, as we have argued 
thoroughly in section 4, this stability criterion actually does not tell us anything about the 
stability of the flux distribution in our case. This is because, for the decaying D-brane, it 
should be noted that the kinetic energy (in addition to the gradient energy) of the tachyon field 
carries the fundamental string charge as well. Consequently, in general, the system does evolve 
even for the total charge held fixed, and this is the physical reason behind the conclusion drawn 
in this work differently from [20, 21, 22]. 

How generic is the conclusion we have drawn? We have analyzed in detail the tachyon rolling 
of the special configuration given in Eq.(2.4), and observed the formation of the string- like thin 
electric flux tube. We also pointed out that modulation of the tachyon field tends to grow by 
using both the boundary-state and the effective field theory descriptions. On the other hand, 
not every inhomogeneous tachyon configuration will eventually evolve to a configuration with 
thin flux tubes like this. In fact, we can easily construct an example with static inhomogeneous 
energy density by turning on the tachyon profile of the form: 

T(x°, x 1 ) = A cosh(:r ) + A' cos^r 1 ). 

This profile also induces an exactly solvable boundary perturbation. Thus, in order to see 
the fate of unstable D-branes with generic tachyon configuration, a more general analysis is 
desirable, which we leave for future study. 

Finally, [13] has proposed an another approach in which the dynamical formation of funda- 
mental string in unstable D-brane is argued within the context of classical effective field theory. 
It would also be interesting to elaborate relation of theirs to our approach. 
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